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Abstract. I summarize recent results Irom radio observations of high redshift star 
forming galaxies, discuss radio continuum emission as a measure of star formation 
rate, and consider future capabilities at cm to IR wavelengths. 

1 Radio surveys to ^tJy sensitivity 

Source counts based on low frequency surveys of the sky with Jy sensitivity 
showed a significant departure from a EucHdean, non-evolving source population, 
indicating, for the first time, cosmic evolution in a source population. The source 
population entailed luminous radio galaxies, with spectral luminosities at 178 
MHz: Pi78 > 10'^^ erg s~^ Hz^^ In these sources the synchrotron radio 

emission is from high energy electrons accelerated in a relativistic jet emanating 



Subsequent observations with sub-mJy sensitivity, starting with the WSRT 
and continuing with the VLA revealed flattening of the source counts be- 
low 5 mJy (Fig. 1). Windhorst et al. hypothesized that this flattening was 
due to a new population of sources, namely star forming galaxies with Pus < 
Iq3i Qj-g g-i Hz~^. The radio emission in these sources is from relativistic elec- 
trons accelerated in supernovae remnant shocks 

There has been a recent revival of deep radio surveys, motivated in large 
part as follow-up to deep optical, infrared, and (sub)mm surveys [|[0,|,|lO). The 
frequency of choice is 1.4 GHz for these deep surveys, allowing for fiJy sensitivity 
with arcsecond resolution and a wide field of view (FWHM = 30'). In the coming 
years the Expanded VLA (EVLA) will push to the sub-/iJy level, while in the 
coming decades the Square Kilometer Array (SKA) will potentially probe nJy 
sources. As pointed out in numerous papers in these proceedings (Bertoldi, Adel- 
berger, Hughes, Sanders), radio observations play an important, complimentary 
role to observations at other wavelengths, in that: 

• They are not plagued by extinction corrections. 

• They provide a rough estimate of source redshifts, in combination with 
(sub)mm observations pT| , p^ . 

• Low order CO transitions redshift into the cm bands, revealing large reser- 
voirs of less dense, cooler gas [ p"3|jl^ . This topic will not be discussed herein. 

• They provide arcsecond astrometry and imaging, thereby avoiding the confu- 
sion problems inherent in deep searches for optical counterparts of (sub)mm 
sources discovered in low resolution single dish bolometer array surveys. 
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For example, at the optical limits of the HDF {Tab < 29) one expects, by 
chance, three faint galaxies within the 6" error circle of a typical SCUBA or 
MAMBO source At 1.4 GHz the source counts between 40 /iJy and 1 mJy 
obey: N(> S1.4) = 2.2 x IQ-^ S^l'"^ arcsec-^ with S1.4 in mJy||. The number 
of spurious Si. 4 > 22/iJy sources within the error circle is only 0.02. 
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Fig. 1. Radio source counts at 1.4 GHz (adopted from Iql). 



Since deep radio images are usually follow-up observations of deep fields at 
other wavebands, there has been very rapid progress in determining the nature of 
the parent galaxies of the /iJy radio source population PJT^ . This short summary 
will focus on recent results on the /iJy radio source population, emphasizing the 
unique information about high redshift star forming galaxies coming from radio 
observations. I will also discuss the radio-to-far infrared (FIR) correlation, re- 
deriving star formation rates based on this correlation, and conclude with a short 
discussion of future instrumentation. 

Richards |^ finds that 75% of the Si. 4 > 40/iJy radio sources are identified 
to / < 24, with a median of / = 22.1. Interestingly, he also finds that 25% 
of the sources are unidentified to / > 25. And perhaps more interestingly, a 
number of groups []6|, pT|Jl^ , |l^ find that (sub)mm observations of these fiJy radio 
sources with faint (or absent) optical counterparts results in a > 40% detection 
rate at mJy levels. This suggests that the optically faint, /zJy radio sources are 
equivalent to the mJy (sub)mm source population. 
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Haarsma et al. [ p^ have used an extensive spectroscopic and photometric 
redshift analysis to determine the redshift distribution of the /xJy radio source 
population. They find that (roughly): 

• 50% of the sources are spirals, or irregular galaxies, at z < 1. 

• 25% are ellipticals, presumably low luminosity AGN, at 0.3 < z < 1.5. 

• 25% are optically faint (or absent), and red. 

They propose that these later sources are likely to be high redshift (1.5 < z < 4), 
dust obscured starbursts. This idea is consistent with the SCUBA and MAMBO 
results discussed above []6|, pj|Jl^ , |l^ . 

The angular size distribution of the /iJy radio source population remains a 
point of debate. WSRT observations of the HDF at 1.4 GHz, 15" resolution, 
to 8/xJy rms detect a number of sources not detected in the VLA survey at 
similar sensitivity but with 1" resolution This would suggest a significant 
population of sources larger than 1". However, confusion is a serious issue at 
this sensitivity level at 15" resolution. Combined MERLIN+VLA observations 
of the HDF suggest that most of the sources have angular sizes between 0.7" 
and 2", with a median of 1.4" (Muxlow et al. in prep). On the other hand, 
VLA observations of the cluster A2125 suggest that most of the sources are 



unresolved, with upper limits of typically 1"|10 . If the sources are indeed 1" in 
size, this presents a significant problem for the SKA, since the sky will become 
naturally confusion limited at the few nJy level, independent of the resolution 
of the instrument [Q, well above the sensitivity of the SKA. 

Richards ^ finds that the mean spectral index for a 1.4 GHz selected sample 
of sources is -0.8, typical of star forming galaxies. Not surprisingly, an 8 GHz 
selected sample shows a flatter mean spectral index of -0.4. 

Barger et al. |l^ find that only 15% of the radio sources with Si. 4 > 25/iJy in 
the Hawaii Deep Field are X-ray sources with 2 to 10 kcV fluxes of: Ix > 1 x 10^^^ 
erg s"^ cm^^. This result is consistent with the idea that the majority of /^Jy 
radio sources are star forming galaxies. 

One problem with deep radio surveys is the limited area covered, such that 
cosmic variance can lead to substantial differences between counts derived for 
different fields. For instance, the counts at the 100 /iJy level in the Phoenix Deep 
Field are a factor two higher than those in the Hubble Deep Field This 
variance can be seen as the increased scatter below 1 mJy in Fig. 1. 



2 The FIR-radio correlation: deriving star formation 
rates from radio observations and the importance of 
inverse Compton losses 

Most recent derivations of star formation rates based on radio observations 
[ p^[|20| , pl| use equs. 21 and 23 in Condon Q to relate the star formation rate 
(SFR) to the 1.4 GHz spectral luminosity. This relationship was derived from 
the supernova rate and the integrated radio luminosity of the Milky Way. 
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An alternate method for deriving the relationship between radio luminosity 
and SFR relies on the tight correlation between FIR luminosity and radio lumi- 
nosity in star forming galaxies Q|. This correlation is remarkable in the small 
scatter observed over at least three orders of magnitude in FIR luminosity. Fur- 
ther, it holds for optical and IR selected samples |Q. This alternate method uses 
spectral synthesis models for star forming galaxies |2^, assuming a fraction, /c, 
of the bolometric luminosity is absorbed by dust and re-emitted in the infrared, 
and then uses the FIR-radio correlation to relate the radio luminosity to the 
FIR luminosity. For a 10* yr continuous starburst, solar abundances, a Salpeter 
IMF from 0.1 to 100 Mq, and fc = 1, this calculation leads to: 

SFR = 5.1 X 10^22 M© yr-i (1) 

where L1.4 is the 1.4 GHz spectral luminosity in W Hz^^. Scaling to the same 
IMF limits, this equation implies a factor 2.5 lower SFR than the equations in 
Condon [|. 

There are a number of uncertainties in both calculations. In the case of Con- 
don's calculation the Galactic supernova rate and radio continuum luminosity 
are both uncertain by at least 50%. For the stellar synthesis model calculation, 
the scatter in the FIR-radio correlation leads to a 50% uncertainty, while chang- 
ing parameters in the starburst model changes the predicted SFR, eg. decreasing 
the starburst age to 10^ years increases the SFR by 50%, and there is the uncer- 
tain fc- It is currently not clear which, if either, calculation is correct. Indeed, the 
relative agreement is encouraging, given the very different metods. And given the 
different conditions in different galaxies (eg. the age of the starburst, the dust 
covering factor, the IMF, ...), it is clear that there will be no globally correct 
relationship, only a statistically most likely one. 

The FIR-radio correlation has a very simple heuristic explanation: both the 
FIR and radio emission relate to massive star formation, with the FIR coming 
from dust heated by the interstellar radiation field, and the relativistic elec- 
trons being accelerated in supernova remnant shocks. However, given the number 
of processes and parameters involved, it remains remarkable, and as yet unex- 
plained, as to why the correlation is so tight [^|j2^. For instance, it has long been 
known that the total radio luminosity of galaxies, both normal disks and nuclear 
starbursts, is an order of magnitude larger than that expected from the sum of 
the supernovae jj,^, although see This requires that the relativistic elec- 
trons be stored in the ISM of galaxies for a timescale, t > Y.'^'^^"" x tswR ~ lO^yr, 

' — Ub.ism ^ ' 

where Ub is the magnetic energy density. 

This leads to the question of the importance of inverse Compton losses off 
the microwave background for relativistic electrons in the ISM of high red- 
shift galaxies. The energy density in the microwave background increases as: 
Umwbg = 4.0 X 10~^^(1 -I- z)'' erg cm^^. This is shown in Fig. 2, along with 
the energy density in the magnetic field in a typical spiral arm . The energy 
density in the magnetic field is larger than that in the microwave background 
to 2 ~ 1. Beyond this redshift, inverse Compton cooling will dominate over syn- 
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chrotron radiation, limiting electron lifetimes, ti.4, and leading to a departure 
from the radio-FIR relation. 
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Fig. 2. The energy density of the microwave background, Umwbg, vs. redshift. Also 
shown is the typical magnetic energy density, Ub, in the arm of a spiral galaxy, and 
Ub in compact nuclear starbursts, and the energy density in the radiation field, U^, in 
starbursts. ti.4 is the radiative lifetime of an electron emitting at 1.4 GHz. 



This is not true, however, for compact nuclear starburst galaxies, ie. systems 
with SFRs > 100 year~^ in regions smaller than a few hundred parsecs 
[ p8| . In these systems the energy density in the magnetic field is thought to 
be almost three orders of magnitude larger than in the disk ||2^ , in which case 
Umwbg only becomes relevant at z > 7. But compact nuclear starbursts raise a 
different, related problem: the energy density in the IR radiation field from the 
starburst itself is larger still than that in the magnetic field. This means that 
inverse Compton cooling should remove the synchrotron emitting electrons on 
fairly short timescales (< 10^ yr), and accentuates the question: why do nuclear 
starbursts follow the radio-FIR correlation? There is a large body of literature 
on this issue Q], but as yet no closure. Hence, we trade one problem, inverse 
Compton losses off the microwave background at high z, for a second, inverse 
Compton losses off the starburst IR radiation field at all z. 

Until proven otherwise, we adopt the radio-FIR correlation as a given, and 
see how it can be used in the study of high z star forming galaxies. 



6 Carilli 



3 Using the radio-FIR correlation to study high z star 
forming galaxies, and what the future holds 

Haarsma et al. have derived the cosmic star formation rate density (SFRD) 
vs. z based on /iJy radio samples. They find a steep rise in the density from 
z = to z = 1, as has been found in optical and IR studies. But they also find 
a systematically higher SFRD at all redshifts by a factor 3 relative to reddening 
corrected optical studies. This suggests that even larger dust corrections are 
needed in optical studies, or that optical studies miss a large population of 
dust-obscured galaxies. However, they use the equations in Condon Q to derive 
SFRs from radio continuum luminosities. If we use equ. 1 above instead, the 
radio derived values agree well with the optical values. 

A second area in which the radio-FIR correlation has been used in the study 
of high z star forming galaxies is as a rcdshift indicator |l|jl|,||,|l|,|2|j3|Jl|l . 
The impetus in this case is the very faint optical counterparts being found for 
most faint (sub)mm sources, thereby precluding follow-up optical spectroscopy 
of a large sample of sources. The radio-FIR method relies on the opposing slopes 
of the synchrotron and thermal dust emission in star forming galaxies. Our most 
recent models use the average observed SED for nearby starbursts to relate 
redshift to the observed spectral index between 350 GHz and 1.4 GHz Figure 
3 shows the model, along with a few sources with known redshifts, including 
sources with AGN and starburst optical spectra. The method is admittedly 
imprecise, especially at high redshift, and there are a few degeneracies in the 
solutions, such as the addition of cold dust or the presence of a radio loud AGN 
but it is the only viable alternative for deriving redshifts for the large 
majority (90%) of the faint (sub)mm sources. The redshift distribution for the 
faint (sub)mm source population as derived using the curves from | p2| can found 
in the contribution by Bertoldi in this volume. 

An important point concerning radio follow-up observations of faint (sub)mm 
sources is the relative sensitivities. Comparing MAMBO images with deep radio 
images shows that 70% (10 of 14) of the > 3.5 mJy (ha) sources at 250 GHz 
have radio counterparts with Si. 4 > 22/.tJy within 3" (Bertoldi, this volume). We 
expect only 0.2 chance coincidences. This result lends confidence to the reality 
of the MAMBO detections. It also implies that the current sensitivity of deep 
VLA fields is well matched to that obtained with mm bolometer arrays, and 
in particular, that there is not a dominant population of very high z sources 
(z > 4), or of low z, cold dust sources. 

Figure 4 shows the expected sensitivity of future cm to IR instruments com- 
pared to the expected flux density of Arp 220 at various redshifts. Overall, the 
next generation instruments are well matched to the expected flux density of 
Arp 220 out to z 2 to 8. Clearly, the ALMA is by far the most sensitive 
telescope relative to the dust spectrum, and will detect low luminosity galax- 
ies to very high z. However, this applies to pointed observations, which have 
a limited field-of-view (FWHM ~ 20"). For surveys of fields larger than about 
15' X 15', the next generation bolometer array cameras operating on large single 
dish telescopes will be competitive with ALMA, as will the EVLA and SIRTF. 
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Fig. 3. The relationship between redshift and observed 1.4-to-350 GHz spectral index 
for an active star forming galaxy, derived from the observed SEDs of 17 low z galaxies 
p^]. The dash line shows the model. Squares are observed values for galaxies with 
starburst spectra. Circles are values for sources with AGN spectra. 

Figure 5 shows the flux density of Arp 220 and M82 vs. z at various frequen- 
cies, relative to the sensitivities of future telescopes. The important point in this 
diagram is the interesting source selection function at (sub)mm vifavelengths: 
typical (sub)mm surveys detect luminous star forming galaxies at essentially all 
redshifts, but they miss completely the low z, low luminosity galaxies. Hence 
(sub)mm surveys result in a very clean, but totally biased, sample of sources. 
Radio observations result in a mixture of low z, low luminosity, and high z, high 
luminosity star forming galaxies, as well as radio loud AGN. 

The closing debate at this workshop contrasted the IR vs. submm vs. optical 
views of high z galaxies and galaxy formation. Each side argued that they detect 
the dominant contribution to cosmic star formation at a specific epoch. Yet, each 
side has a very specific galaxy selection function, and the overlap between the 
populations apparently is small, ~ 10%. Indeed, one might argue that /iJy radio 
surveys are the least biased means of detecting all the source populations. But 
therein lies the fundamental problem: how to differentiate the source populations 
on a deep radio image? Overall, it is clear that all sides are currently seeing 
only limited, and perhaps orthogonal, aspects of galaxy formation. In order to 
address the general question of galaxy formation, or at least the formation of 
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the stellar content of galaxies, requires wide field surveys using the EVLA, the 
next generation bolometers arrays on the LMT and GBT, and SIRTF, with very 
deep follow-up studies of selected samples of sources using ALMA and NGST. 

The National Radio Astronomy Observatory is operated by Associated Univ. 
Inc., under contract with the National Science Foundation. I would like to thank 
A. Hopkins for allowing me to reproduce Fig. 1, and B. Poggianti, N. Miller, A. 
Blain, M.Yun, and F. Owen for useful discussions. 
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Fig. 4. The observed spectrum of Arp 220 at various redshifts, compared to the sen- 
sitivities for various existing and future telescopes. 
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Fig. 5. The flux densities of Arp 220 and M82 vs. z at various observing frequencies. 
Also shown are the sensitivities of existing and future telescopes. 
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